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The three-dimensional study of the ring vortex solitons is conducted for both attractive and 
repulsive BECs subject to harmonic potential confinement. A family of stationary ring vortex 
solitons, which is defined by the radial excitation number and the winding number of the intrinsic 
vorticity, are obtained numerically for a given atomic interaction strength. We find that stabilities 
of the ground and radially excited states of the ring vortex soliton are dependent on the winding 
number differently. The ground state of the ring vortex soliton with the large winding number is 
unstable dynamically against random perturbation. The radially excited state of the ring vortex 
soliton with large winding number corresponds to the increased collapse threshold, and therefore can 
be made stable for sufficiently small atomic interaction strengths. The ground and radially excited 
states also demonstrate different dynamical evolutions under large atomic interaction strengths. 
The former exhibits simultaneously symmetrical splitting in the transverse plane, while the latter 
displays periodically expand-merge cycles in the longitudinal direction. 

PACS numbers: 03.75.Lm,05.45.Yv 



I. INTRODUCTION 

The realization of Bose-Einstein condensates (BECs) 
in dilute quantum gases has drawn a great deal of in- 
terest in vortices [U E] • Recently, considerable efforts are 
aimed at the prediction of settings supporting stable mul- 
tidimensional solitons with intrinsic vorticity. One such 
structure is the ring vortex soliton, a soliton loops back 
on itself to form a ring and the phase of the wave function 
winds through an integer multiple of 2tt radians around 
the vortex line. Such ring vortex soliton can be defined 
by the quantum number set (n, 5), where n is the radial 
excitation number, equal to the number of the rings, and 
S is the winding number of the intrinsic vorticity. One 
example of the ring vortex soliton, being called "ring- 
profile solitary waves" , has been studied in the context of 
optics H3IU- Another example of the ring vortex soliton, 
called "vortex tori" , has also been predicted in three di- 
mensions (3D) cubic-quintic Ginzburg-Landau equation 
[5j |6] . The ring dark soliton is first introduced in the re- 
pulsive BECs [7 j. In the recent binary BEC experiment, 
the ring bright solitons are also observed [8j [9] , which is 
studied further in the theoretical works [TU] . 

The most simple structure of the ring vortex soliton is 
of one ring in the radial direction, i.e., n=l, which can be 
regarded as the radial ground state, and here we call it as 
single ring vortex soliton (SRVS). The bright vortex soli- 
ton created in the attractive and dipolar BECs [TTHT6] 
can be regarded as one typical example of SRVS. As we 
are known, the atomic interaction in BECs can be con- 
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trolled by the Feshbach resonance technique [T7HT9] to 
change the strength and sign of the interaction. There- 
fore, great interests are motivated to investigate the prop- 
erties of the vortex soliton under the nonlinearity of the 
atomic interaction. It is predicted that SRVS is stable 
under the 3D harmonic trapping potential for sufficient 
small interaction strength [TTJ [T2l [20] , while for large in- 
teraction strength the dynamical instability phenomena, 
such as split-merge cycles of particles [TTJ [T4J EI] and in- 
tertwining of doubly quantized vortex [22H24] , occur in 
the unstable regime. 

An infinite sequence of radially excited stationary 
states of the ring vortex have been predicted to Gross- 
Pitaevskii equation (GPE) in two dimensions (2D) and 
sphere shell in three dimensions (3D) [25j[26]. The ra- 
dial excited state is of multiple concentric density-wave 
rings, i.e., n > 2, and here we call it as the multiple ring 
vortex soliton (MRVS). The one dimensional (ID) and 
2D solutions of the MRVS are also proposed for a spe- 
cific spatially- modulated nonlinearity [27-29 . In experi- 
ments, the ringlike excitations have been observed in the 
hyperfine states \F = 1, rrif = —1) and \F = 2, rrif — +1) 
of 87 Rb under the rotating cylindrically magnetic trap 
[5]. In | ± 1) spinor 8T Rb BECs, more ringlike excita- 
tion modes have also been yielded under cylindrically 
magnetic confinement [9 . In theoretical simulations, a 
family of 3D gap solitons of the multiple rings has been 
reproduced, as supported by ID optical lattices [30] . 

In addition to the intrinsic vorticity, the MRVS has 
an additional excitation freedom, corresponding to an 
infinite number of nodes of the wave function in the ra- 
dial direction. The early linear stability analytical results 
[25j[26] have shown that, despite of the radial excitation, 
the 2D (or 3D in the sphere shell) MRVS can be stable 
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in the harmonic confinement below a threshold of the 
interaction strength. For large interaction strength, the 
MRVS states are unstable to collapse. When the instabil- 
ity times is much larger than the time scale 2tt /uj, uj is the 
harmonic trap frequencies, it is said to be experimentally 
stable for small interaction strength [25, 26 . However, in 
a recent 3D study, it is found that the MRVS with n=2 
is unstable against quadrupole perturbations [30]. This 
implies that the radial excitation of MRVS could display 
more rich stability properties if subjected to the 3D ge- 
ometry. This motivates us to investigate the stability 
of the MRVS in the 3D space. Besides, considering the 
cylindrical magnetic confinement in the related experi- 
ments [H |9] , the full 3D equation is necessary to study 
the ring vortex solitons. 

In this paper, the 3D ring vortex solitons are con- 
structed in the framework of GPE at the cylindrical co- 
ordinate for both attractive and repulsive BECs. The 
stabilities of SRVS and MRVS depend on S differently. 
The SRVS with S=l has better stability, i.e., large g c , a 
threshold of the atomic interaction strength below which 
the ring vortex state is robustly stable against random 
perturbation, while the SRVS with S > 2 corresponds to 
the greatly decreased g c and therefore poor stability. On 
the contrary, the MRVS with large S has better stability 
especially when the radial excitation is high. No or very 
small stable regime is found for the MRVS with S=l. 
The differences between the ground and radially excited 
states of the ring vortex soliton are further demonstrated 
in the dynamically unstable evolution. The SRVS shows 
the simultaneous symmetrical splitting in the transverse 
plane, while the MRVS shows the periodical expand- 
merge cycles in the longitudinal direction before collapse. 
Such dynamical instability undergoes at the timescale of 
about several seconds, being longer than 2tt/uj, proving 
their good experimental stability. 

The paper is outlined as follows. We first give a fam- 
ily of ring vortex solitions by using Newton continuation 
method in Sec. II. Next, the stability properties of the 
ring vortex solitons is analyzed by using linear stability 
analysis in Sec. Ill, which is followed by the direct nu- 
merical simulations of the perturbed ring vortex solitons 
in Sec. IV. Finally, we conclude the main results of the 
work in Sec. V. 



II. THEORETICAL MODEL AND STATIONARY 
SOLUTIONS 

We consider the BECs in an external harmonic trap- 
ping potential V(r, z)=m(uj 2 r 2 +uj 2 z 2 )/2, where r 2 =x 2j r 
y 2 , mis the atom mass, uo r and uo z are the radial and axial 
trapping frequencies. The wave function ip of the BECs 
satisfies the dimensionless GPE 
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where 6 is the azimuthal angle, g=4irNa s /ao is the 
interaction strength, as determined by the total num- 
ber of particles N in the condensate, the s-scattering 
wave length a s , and the harmonic oscillator length 
ci(j=\/h/(jjm : Q is the rotation angular frequency. The di- 
mensionless V(r, z) = \(^ 2 r 2 +j 2 z 2 ) with7 r =^, 7^ = ^, 
and cj=min{coV, uj z }. Eq.(l) is obtained by rescaling the 
length by ao, the time by Co> -1 , and the energy by hw. 
Assume the wave function is 
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where S is azimuthal quantum number, i.e., the intrinsic 
vorticity and \i is the chemical potential. The function 
<p(r : z) satisfies the equation 

(3) 

with the boundary conditions of lim <fr(r, z) = and 

r— >-0 

lim 0(r, z) = 0. The linear-limit solution of Eq. (3) is 
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written as the scaled linear combinations of products of 
the wave functions of the harmonic oscillator [31], 

(5+l)/2 1/4 

<f>(r,z) = lr * r s L s n ( lr r 2 )e-^W), (4) 

where C% = nf=i( n + i)> ^ s ^ ne generalized-Laguerre 
polynomials, and n is the quantum number in the r di- 
rection (the quantum number in the z direction has set 
to be zero). The chemical potential corresponding to Eq. 
(3) and Eq. (4) is written as 

= 7r[2(n - 1) + S + 1] - SQ + )p z . (5) 

The stationary soliton solutions can be obtained as nu- 
merical solutions of Eq. (3) using the Newton contin- 
uation method [32] with fi n s input, starting with the 
linear-limit solutions, under the constraint that the nor- 
malization of the wave function N and the energy E are 
conserved with numerical iteration time, where 

7V(V0= / |^(R,t)| 2 ^R = ^o) = i,t>o, (6) 

and the energy 
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Shown in Fig. 1 are various ring vortex solitons ob- 
tained in the case of the attractive interaction under dif- 
ferent shaped harmonic potentials. For the isotropic har- 
monic potential of j r =j z =l, g=-0. 00837V refers to a T Li 
system both with coV=co^=207rHz, for the cigar-shaped 
harmonic potential of 7 r =1.44 and 7^=1, g=-0. 01177V 
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FIG. 1: The 3D ring vortex solitons with different quantum- 
number sets (71,5) found numerically for attractive BEC, as 
visualized by isosurface of dimensionless density of the wave 
function in the isotropic (a)-(h), pancake (i)-(l), and cigar 
(m)-(p) harmonic traps, respectively. 



refers to 7 Li with u; r =178Hz and co^=123Hz [33], and for 
the pancaked harmonic potential of 7 r =l and 7^=20, g=- 
0.0083iV refers to 7 Li with cj r =207rHz and c^=400ttHz 
[34] . In the case of isotropic harmonic potential, we 
also study the ring vortex soltons for the repulsive in- 
teractions, where #=0.01887V refers to 87 Rb BECs with 
coV=co^=207rHz. In this work, Q=0.7 is fixed. The vor- 
tex solitons obtained for the repulsive interactions (not 
shown here for the sake of brevity) resemble those shown 
in Fig. 1 for the attractive interaction. 

The ring vortex solitons are further illustrated accord- 
ing to the radial profile of the wave function. Fig. 2(a) 
shows that the influence of the quantum number n on 
the ring vortex solitons, i.e., the radially excited level is 
increasing with n, corresponding to more nodes of the 
wave function in the radial direction. Fig. 2(b) shows 
that the influence of the quantum number S on the ring 
vortex solitons. The ring vortex solitons with the larger 
S demonstrates the deeper depletion of density at core 
area (denoted by dash arrows) and the lower density peak 
(denoted by solid arrows). Fig. 2(c) shows the influence 
of the sign of the interaction nonlinearity on the shape of 
the ring vortex solitons for a given interaction strength. 
The density of the wave function at the first peak (de- 
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FIG. 2: The radial density profile of the ring vortex soli- 
tons characterized by the quantum numbers (n, S) for the 
attractive (g < 0) or repulsive (g > 0) interactions under 
the isotropic harmonic potential in (a)-(c) and also the cigar- 
shaped and pancake-shaped harmonic potentials in (d). 



noted by the arrow in Fig. 2(c)) in the case of the at- 
tractive interaction is higher a little than that in the case 
of the repulsive interaction. Such difference is the most 
prominent in the case of (n, 5) =(1,1), but is weak in the 
cases of large n or large S. Fig. 2(d) illustrates variation 
of the radial profile of the ring vortex solitons under the 
different shaped harmonic potentials. 

Shown in Fig. 3 is the /jl(N) curves of the ring vor- 
tex solitons. For a given TV, as compared with the SRVS 
family with n=l, the MRVS families with n=2, 3, 4, ... 
are of increasing \i. Within a SRVS or MRVS family, \i 
is increasing with S. Each data shown in Fig. 3 corre- 
sponds to a convergent stationary state, which is further 
investigated for the stability and dynamical evolution. 



III. STABILITY ANALYSIS 

The stability of the ring vortex solitons is analyzed 
using the linear stability analysis. The wave function 
that deviates slightly from the stationary solutions is con- 
structed as 

^ = [0(r, z) + ue iEt + w""e- iEH }e iSe -^\ (8) 

where \u\, \w\ <C 1 are eigenmodes. Substituting Eq. (8) 
into Eq. (1) gives 
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Sfl + /i, E is the eigenvalue related to u and w, which is 
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FIG. 3: (Color online) The /x(iV) curves of the ring vortex soli- 
tons of n=l ,2, 3, 4 with 5=1 (square), 5=2 (circle), and 5=3 
(triangular), as calculated with j r =j z = l and g=-0.0083iV, 
and Q=0.7, referring 7 Li condensate under the isotropic har- 
monic potential with cj r =c<;; Z =207rHz. 



obtained by diagonalization of Eq. (9) under the bound- 
ary conditions demanding that u(r,z),w(r,z) — >• at 
r,\z\ — >• oo and r, \z\ — >• 0. Applying Eq. (9), we analyze 
the ring vortex solitons of different radial excited states, 
not only the SRVS, i.e., n=l, but also the MRVS, from 
the first radially excited state, i.e., n=2, to the higher 
radially excited states. 

Fig. 4(a) shows the real and imaginary parts of 
eigenvalues of the SRVS obtained for a given interac- 
tion strength in an isotropic harmonic potential. When 
5=1, all eigenvalues are real, indicating the correspond- 
ing state is stable. When 5 >2, complex eigenvalues 
emerge, indicating the corresponding states are dynami- 
cally unstable. The inverse of the imaginary part of the 
complex eigenvalue, Im E, gives the time scale of such a 
dynamical instability. As shown in Figs. 4 (a) and 4(b), 
it is found that Im E decreases with 5, which indicates 
that the ring vortex soliton with the larger intrinsic vor- 
ticity is of the longer lifetime before collapse, i.e., the 
better experimental stability. 

Moreover, Figs. 4(c) and 4(d) show the g dependence 
of Im E obtained for the SRVS with different 5. The 
eigenvalues become complex above the critical interac- 
tion strength g c , showing the onset of the dynamical in- 
stability. The g c value is dependent on 5. For example, 
our results show that g c is around -11 when 5=1 but 
decreases to -1.5 and -2 when 5=5 and 20 in the case 
of the attractive interaction. The corresponding value 
of g c becomes a little larger in the case of the repulsive 
interaction, being around 12, 3, 1.5 for the SRVS with 
5=1, 5, and 20, respectively. The case of the SRVS, 
i.e., the general ring vortex soliton without radial exci- 
tation, has been studied previously (TTJ [I4j |20] . It was 



0.3 
0.2 

lu - 1 
E 

_ 0.0 
-0.1 
-0.2 
-0.3 



I 1 

(a) n-i , 


i ■ i 1 

U— - I I 

a O 

A 

O 


i 

□ 

o 

A 

v 
O 


S=1 • 
S=2 - 
S=5 
S=10- 
S=20 










€ 


) ! 

V O 






' (1,1) t 


=342 









,90=3^.^-90=12/ 
V 



12 -16 3 6 



0.10 - 





i j • i / 


0.10 




1 ' 1 f 1 


(e) n=2, g<0 


/ 

J / 

/ * 

/ / 

/ / 


0.08 

W0.06 
E 


• (f) n=2, g>0 

/ / 

/ .lL. 


/ 9 

9c= " 10 //." 


A' 

• A * 

• •▲• • S=20 


0.04 
0.02 
0.00 

( 

0.08 


g c =1.£ 


■" . ' ' 

if'' -B-S=1 
P -•■ S=5 ■ 

■•A--S=20 


-4 
> 1 < — 


-8 -12 

g 

— 7—1 > r»- 


) 


4 8 12 

g 

1 ' 1 / ' 1 




FIG. 4: (Color online) The imaginary part in energy spec- 
trums (Im E) of the stability analysis, (a)-(b): Im E .vs. 5 
obtained for SRVS under a given attractive interaction g=-ll, 
inset in (a) displays a robust SRVS wit 5=1 at dimensionless 
time t=342 against noise perturbation, and Im E .vs. g ob- 
tained for SRVS [(c)-(d)] and MRVS of n=2 [(e) and (h)] and 
n=4 [(g) and (h)j with different 5. 



found that for sufficiently weak interaction strength the 
SRVS with 5=1 is stable, while for 5 >2 the SRVS is 
unstable to quadrupole oscillations. For the instability 
times being much longer than the time scale 2tt/uj of the 
BECs, the SRVS with 5 >2 is said to be experimentally 
stable for small interaction strength. Our results shown 
in Figs. 4(a)-(c) are consistent with the previous study 
for the general ring soliton without the radial excitation 
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Figs. 4(e) and 4(f) show the curve of Im E .vs. g ob- 
tained for the MRVS of n=2, i.e., the first excited state. 
The stability of the ring vortex soliton of the first ra- 
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FIG. 5: (Color online) The dynamical evolution of the MRVS 
in unstable and stable regions in Fig. 4 (e) and (f), (a) is dy- 
namically unstable and (b)-(e) are dynamically stable states, 
which are consistent with linear stability analysis. Referring 
a 7 Li system confined by uJr—uJz—2^ Hz, the size is about 
120/xm and the real time £'=(). 0159£ second, t=0 corresponds 
to the initial state of a 1% random perturbation in density. 



dial excited state has been studied based on the 2D solu- 
tions [25 , which shows that the first radial excited state 
with 5=1 is stable for the sufficiently small interaction 
strength below a critical value. In our 3D study, for 5=1, 
no stable regime is found for the first excited state. This 
result means that the 3D MRVS of the first excited state 
with 5=1 could be different from its counterpart in 2D. 
The critical dimensionality for the GPE is 2D. Further, 
we find that, for large 5, the 3D MRVS can be made 
stable. For example, we gives two Im E .vs. g curves ob- 
tained for MRVS with 5=5, and 5=20 in Figs. 4(e) and 
(f), where the stable regime, though being very small, is 
found up to g c =-1.0 in the cases of the attractive interac- 
tion and g c =1.5 in the case of the repulsive interactions. 
Such results indicate that the stability of the MRVS is 
dependent on the intrinsic vorticity. The MRVS with 
large 5 is expected to stable. 

Shown in Figs. 4(g) and 4(h) are results for the MRVS 
of n=4. The 5 dependence of the stability becomes more 
prominent when the MRVS is of the higher radial excita- 
tion. The stable region, though being very small, below 
g c =-0.2 and # c =0.4, are found for the n=4 MRVS with 
5=1 under the attractive and repulsive interactions. The 
stable region increases for large 5. When 5=5, the stable 
regions are up to g c =-2 and g c =3 in the cases of the at- 
tractive and repulsive interactions respectively. Further, 
when 5=20, the stable regions increase up to g c =-^ and 
g c =5 in the cases of the attractive and repulsive interac- 
tions respectively. Through comparing two cases of n=2 
[Figs. 4(e)- 4(f)] and n=4 [Figs. 4(g)-4(h)], we can find 
that the intrinsic vorticity tends to improve the stabil- 
ity of the MRVS, and such effect turns more significant 
when the radial excitation of MRVS is high. The MRVS 
of large 5 has the prominent stable region. 

The linear stability analytic results for the 5 depen- 
dence of the MRVS stability is further confirmed by the 
directed simulation based on Eq. (1). Corresponding to 
Fig. 4(e), Fig. 5(a) and (b) proves that the radially first 
excited state of MRVS of n=2 is unstable against random 
perturbation when 5=1, but turns to stable when 5 is 
large, such as 5=5, under a small g. Next, correspond- 
ing to Fig. 4(g), Fig. 5(c)-(e) prove the stable MRVS of 
n=4 with different 5 under the given interaction strength 
below the increasingly large g c . 



IV. DYNAMICAL EVOLUTION 

The ring vortex solitons are unstable against the dy- 
namical instability when the interaction strength is larger 
than g c . Such dynamical instability is simulated through 
numerically integrating Eq.(l) using the time- splitting- 
spectral technique [35 . In our simulation, the stationary 
solutions of the ring vortex solitons are used as initial 
states after adding a random perturbation of a 1% rela- 
tive amplitude of density. 

We first simulate the dynamical evolution of the SRVS. 
As shown in Fig. 6, the SRVS evolves into polygonal 
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FIG. 6: (Color online) Noise-induced splitting of the ring vor- 
tex solitons of different S at the dimensionless time £, and 
the corresponding density and phase (the third and forth 
columns) plots in the x-y plane. To translate the results into 
the experiment-related units, we assume a 7 Li condensate con- 
taining about 1500 atoms in a isotropic harmonic potential 
with cj r =uj z =207r Hz. The radius of the ring vortex solitons 
here are about 120/im, decaying at the real time £=0.0159t 
second. 



rings, then splits into small fragments, and finally col- 
lapses. The multiple symmetry is revealed during such 
evolution. As S increases, we obtain the 2-, 3-, 4-, 5-, 6-, 
7-, and 8- fold symmetry, as shown in Figs. 6 and 7. As 
seen in the phase profile of the wave function, the multi- 
ply quantized vorticity splits into many singly quantized 
vortices. Such singly quantized vortex cannot be seen 
in the density plot, and hence they are called as hidden 
(or ghost) vortices [36j[37]. The singly quantized vortices 
are self-organized into the regular lattice, contributing to 
the symmetry in splitting. Such symmetry in instability 
can be regarded as simultaneous one, because the ran- 
dom initial perturbation we add is of no azimuthal mode. 
The multiple symmetrical splitting phenomena have been 
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FIG. 7: (Color online) Noise- induced splitting of giant vor- 
tices. The splitting numbers as a function of different S for 
the attractive interaction are shown. 



demonstrated for giant vortices [38-41 j. The symmetry 
in splitting and the number of the fragments are consid- 
ered to be equal to the added angular-momentum quan- 
tum number of the Bogoliubov excitation mode respon- 
sible for the splitting [38, 41 . Ones find a total of three 
types of splitting patterns, 2-fold, 3-fold, and 4-fold sym- 
metries, when the initial perturbation is random noise in 
the density of the wave function [40j|4T]. Figs. 6 and 7 
indicate that the symmetry in splitting of the ring vortex 
soliton can have more types. As matter of fact, the num- 
ber of symmetry can be manipulated by the interaction 
strength g. As shown in Fig. 8, the number in symmetry 
(or the number of the fragments) decreases with g in the 
case of the repulsive interaction and increases with g in 
the case of the attractive interaction. 

Second, we simulate the dynamical evolution of the 
MRVS. As shown in Fig. 9, the MRVS demonstrates 
expand- merge cycles, then collapses. The similar phe- 
nomenon, expand-shrink cycles occurring in the trans- 
verse direction, was predicted for the dark SRVS [TTJ [I4j 
l2Tj . As compared with the dark SRVS, a principal dif- 
ference is that the expand-merge cycles of the 3D MRVS 
occurs in the longitudinal direction. The condensate ex- 
pands along the axial direction. The outer rings expand 
faster than the internal rings. At £=54, the outmost ring 
splits into two pieces, and the internal rings expand but 
remain united. In contrast, the inmost ring shrinks, as 
seen in the density profile in the longitudinal plane shown 
in Fig. 9, middle row. The expanded condensates and 
pieces subsequently unite to restore the original shape at 
£=162. Such cycles repeats two times, then collapse fi- 
nally. The expand-merge cycles reflect the expansion in 
the axial direction and the oscillation in the radial direc- 
tion of the density of the wave function. Besides, during 
the dynamical evolution of the MRVS, as seen in the 
phase profile of the wave function in the transverse plane 
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shown in Fig. 9, lower row, the multiply quantized vor- 
ticity is maintained well during the expand-shrink cycles. 
No vortex splitting is observed until the final collapse, 
which indicates that for the MRVS, the instability of the 
multiple radial waves prevails over the azimuthal insta- 
bility. The similar expand-merge evolutions are demon- 
strated in the case of the repulsive interaction, where the 
cycles repeat three times before collapse, indicating the 
longer lifetime. 

We also simulate the dynamical instability under the 
pancake-shaped and cigar-shaped harmonic potentials. 
The expand-merge cycle occurring for the isotropic har- 
monic potential is suppressed under the pancake confine- 
ment potential. Under the cigar-shaped harmonic po- 
tential, being elongated in the axial direction, the MRVS 
does not demonstrate expand-merge cycles either. Be- 
sides, another interesting dynamics phenomenon, inter- 
twining of vortices in the cigar-shaped BECs [22-24 , is 
also not observed in our simulations. 



FIG. 8: (Color online) Noise- induced splitting of the 3D vor- 
tex solitons of large vorticity in (a)-(c) with S = 25 for the re- 
pulsive interaction and (d)-(f) with 5=45 for the attractive in- 
teraction. To translate the results into the experiment-related 
units, we assume the 87 Rb condensate containing about 640 
(a), 4360 (b), and 8000 (c) atoms and the 7 Li condensate con- 
taining 1325 (d), 4820 (e), and 9638 (f) atoms in a isotropic 
harmonic potential with uj r —(^ z —^07i Hz. The radius of the 
ring vortex solitons here are about 31/xm for (a)-(c), and 
120/im, being snapshots at the real time t' = 0.0159t second. 
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FIG. 9: (Color online) Noise-induced splitting of a vortex 
soliton of (n,S) = (4,3) at the dimensionless time t and the 
corresponding density plot in the x-z plane (middle row) and 
phase plot in the x-y plane (lower row). To translate the 
results into the experiment-related units, we assume a 7 Li 
condensate containing about 1500 atoms in a isotropic har- 
monic potential with uj r =w z =207r Hz. The radius of the ring 
vortex solitons here are about 83.68/im, decaying at the real 
time £ =0.0159£ second. 



V. DISCUSSION AND CONCLUSIONS 

Consider the cylindrical magnetic confinement, we use 
fully 3D equation to investigated the stationary state, 
stability, and dynamical evolutions of the ring vortex soli- 
tons in the attractive and repulsive BECs. A family of 
stationary ring vortex solitons are obtained numerically 
by Newton continuation method. The stability proper- 
ties of the ring vortex solitons are predicted for a given 
interaction strength by using the linear stability analysis 
and confirmed further by the direct simulation. The sta- 
bilities of the SRVS and MRVS depend on 5 differently. 
The SRVS with 5=1 corresponds to a large g c , a thresh- 
old below which the solution is stable against random 
perturbation, while the SRVS with 5 > 2 corresponds 
to the greatly decreased g c . The 5 dependence of the 
stability of the MRVS is on the contrary. The prominent 
stable regimes are found for the MRVS of n=4 with large 
5, such as 5=5 and 20, while very small stable regimes 
are obtained for the MRVS of n=4 with 5=1, and even 
no stable regime is found for MRVS of n=2 with 5=1. 

Therefore, we can expect the robust dynamical stabil- 
ity against random perturbation for the SRVS with 5=1 
and the MRVS with large 5, when the atomic interaction 
strength is less than g c . For the radial ground state of the 
ring vortex solition, the value of g c are around -11 and 12, 
as predicted for the SRVS with 5=1 under the attractive 
and repulsive interactions. For the radially excited state, 
the values of g c are around -4 and 5, as predicted for the 
MRVS with 5=20 under the attractive and repulsive in- 
teractions. In our study, we consider the cylindrical mag- 
netic trap with the harmonic frequency CoV=cj 2 =207tHz, 
which corresponds to a critical atom number of the sta- 
ble ring vortex soliton being about 600 ~ 1300 in the 
case of the radially ground state (where, we refers to the 
87 Rb system with a s =5lA and ra=1.44x 10 _25 kg and the 
7 Li system with a s =-79.35A and ra=1.1702xl0" 26 kg, re- 
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spectively) and several hundreds in the case of the radi- 
ally excited state. To decrease the magnetic harmonic 
frequency down to several Hz, one can expect that the 
critical atom number increases to 10 3 ~ 10 4 . Our re- 
sults suggest the possible conditions that the ring vortex 
soliton can be make stable against random perturbation. 

On the other hand, our study also suggest some most 
unstable states of the ring vortex solitons, for example, 
the SRVS with S > 2 and the MRVS with 5=1. Our di- 
rect simulations show that the lifetime of the dynamical 
evolution can be several seconds before the final collapse, 
for example the SRVS shown in Fig. 6 and the MRVS 
shown in Fig. 9 contain about 10 3 atoms can survive 
up to 1.5-8 seconds if referring to a magnetic trapping 
potential with frequency uj=20tt. Such timescales are 
far larger than the time scale 2tt/uj [25, 26 j. Therefore, 
the dynamical evolution results further suggest that the 
ring vortex solitons could be observed in experiments if 
subjected an appropriate interaction strength. The dif- 
ferences between the ground and radially excited states 
of the ring vortex soliton are also reflected in the different 
evolution ways before collapse. In the dynamically unsta- 



ble regime subjected to the large atomic interaction, the 
SRVS demonstrates the simultaneous symmetrical split- 
ting in the transverse plane, while the MRVS exhibits 
the periodical expand-merge cycles in the longitudinal 
direction. 

Our results suggest possibilities for creation and obser- 
vation of robust 3D ring vertex solitons in the cylindri- 
cal geometry under the magnetic harmonic confinement. 
The present work is conducted under the rotating GPE 
frame, which facilitates the ongoing study that considers 
the possible dependence of the stability and dynamical 
evolution of the ring vortex soliton on the rotation angu- 
lar frequency. 
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